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Abstract 

We develop a new hadronic model for GEANT4 that is specialized for the disintegration of the 
deuteron by photons, d'j —>• np. For the description of two-nucleon interactions, we employ a 
pionless effective field theory with dibaryon fields (dEFT). We apply the new model of GEANT4 
(G4dEFT) to the calculations of the total and the differential cross sections in dy —)• np and 
compare the results with empirical data. As an application of the new model, we calculate the 
neutron yield from the 7 -|- GD 2 process. G4dEFT predicts peaks for the neutron yield, but the 
existing model of GEANT4 does not show such behavior. 
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I. INTRODUCTION 


Effective field theory (EFT) has become a popular method to study hadronic reactions 
with and without external probes at low energies. In the low-energy region where the scale 
of momenta is much smaller than the mass of pions, treating the exchange of pions, as well 
as heavy degrees, in terms of effective contact interactions may be reasonable. Theories 
thus constructed are called pionless effective held theory SHa. In this work, we employ a 
modihed version of pionless theory, so-called pionless theory, with dibaryon helds (dEFT) 
[S]. In the dEFT, a perturbative expansion is directly manifested in physical observables 
such as cross-sections and their analytic forms are available. 

Particle transport codes such as GEANT4 (GEometry ANd Tracking) [9l HO] are com¬ 
monly utilized in the setup of plans for efficient experiments and in the analysis of results. 
Doing so, we can accurately estimate observables for the reactions under consideration. How¬ 
ever, one of the authors (JWS) recently reported that built-in hadronic models of GEANT4 
(vlO.O) fail to describe peaks that are produced through the ®Be(p,n)®B reaction in the low- 
energy region m- GEANT4 was originally built for high-energy physics. In some cases, 
it fails to describe low-energy phenomena correctly as was shown in the ®Be(p,n)®B reac¬ 
tion. In Ref. 11, the authors developed a new data-based charge-exchange model for the 
®Be(p,n)®B reaction, and were able to achieve good agreement with experimental data. 

Inspired by the observation in Ref. 11, we extend the application of GEANT4 models to 
low-energy processes in the two-nucleon systems. We simulate the dj —)■ np reaction by using 
GEANT4 code and hnd that an existing model of GEANT4 (vlO.l) gives a null result for the 
total cross-section of the dy —)■ np reaction at energies below the pion threshold. In order to 
obtain a better description at energies below the pion threshold, we construct a dEFT-based 
hadronic model of GEANT4 for the dy —)■ np reaction. To check the validity of our model, 
we simulate the total and the differential cross-sections of the dy —)■ np process by using 
GEANT4 with dEFT (G4dEFT) at various energies and angles, and we compare the results 
with available experimental data. The result proves the usefulness of the combination of 
dEFT with GEANT4. As an application of the model, we calculate the neutron yields from 
the y -|- GD 2 reaction. We hnd that the G4dEFT model predicts results critically different 
from those predicted by using existing GEANT4 models. 

We organize the paper as follows; In Section II, we briehy present analytic formulae for 
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the differential cross-section of the dEFT, and the description of GEANT4 models relevant 
to this work. In Section III, we present the results for the cross sections in dj —)■ np with 
G4dEFT and compare them with experimental data. We propose experiments from which 
we can test the predictive power of G4dEFT. We summarize the work in Section IV. 


II. METHODS 


Large scattering lengths signal the presence of weakly-bound or almost-bound states. 
Through the inclusion of dibaryon fields that represent a weakly-bound state {^Si channel) or 
an almost-bound state (^S'o channel) of two nucleons, this formulation simplifies calculations 
compared to the pionless theories without dibaryon fields [T2IIT4] . We showed that the 
dEFT up to next-to-leading order (NLO) could be applied to diverse two-nucleon systems 
successfully at low energies [T4I12U] . 

The differential cross-section for the dj —)■ np reaction can be written as 
da a pEi 


dQ 247r k 




( 1 ) 


spin 


where a is the fine structure constant, k is the energy of an incoming photon, and Ei = 
rn^j -f p^ is the energy of an outgoing nucleon in the center-of-mass (c.m.) frame. The 
spin-averaged square of the amplitude can be written as 


S = 1Q{\Xms\^ + \Ymv\^) + ^{\Xmv\^ + \Yms\^) 

spin 

+ l2[l-{k-pf]{\XE? + \YE?), (2) 


where S' is a symmetry factor for the spin average, S = 2. The matrix elements Xmv-, Xms, 
Xe, Ymv, Yms, Ye, and the definitions of k and p can be found in Ref. 20. 

GEANT4 (GEometry ANd Tracking) is a simulation tool kit written in the G-I--I- lan¬ 
guage, which allows microscopic Monte Garlo simulations of the propagation of particles in¬ 
teracting with materials. It is being widely and successfully used in many different scientific 
fields, such as neutron shielding studies [2I1I22], medical physics [23H2S], accelerator-based 
radiation studies [II1126H2H], environment radiation detection studies [29H3lj, etc. 

There are many GEANT4 physics models and cross-sections for both hadronic 
and electromagnetic (EM) interactions. For photonuclear interactions, however, the 
only GEANT4 cross-section and hadronic model relevant to this work are the 
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FIG. 1: (Color online) Total cross-section of the d'y — )• np reaction with respect to the incident 
photon energy in the laboratory frame {E^). Symbols denote the experimental data |34H42] taken 
from the EXFOR database [l5]. Dot-dashed and dashed lines represent the nuclear data extracted 
from ENDF/B-VII.l [33] and JENDL/PD-2004 [31|, respectively. Black and red solid lines are the 
results obtained from G4 and G4dEET, respectively. 

“G4PhotoNuclearCrossSection” class [52] and the “G4CascadeInterface” [55] class, respec¬ 
tively. The G4PhotoNuclearGrossSection gives the total inelastic cross-sections for pho- 
tonnclear interactions. A Bertini-style cascade model is used for G4GascasdeInterface to 
calculate the hnal states of photonuclear reactions. The Bertini model is available for the 
incident photon energies < 3.5 GeV in GEANT4 (vlO.l), but the model has been 

tested mostly in the energy range 60 MeV ~ 3 GeV [55]. GEANT4 uses them by using the 
“G4EmExtraPhysics” GEANT4 Physics Gonstructor. 


III. RESULTS 

A. Total Cross-Section 

To check the validity of GEANT4 for the dy —)■ np reaction, we hrst calculate the total 
cross-section of the reaction by using G4EmExtraPhysics (G4). In Fig. we show the result 
of G4 as a black solid line, which is consistently zero in the energy range < 70 MeV. 
Evidently, the simulation using GEANT for the few-nucleon systems at low energies must 
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FIG. 2: Ratio of the calculated total cross-section for the d'y —)• np reaction obtained by using 
G4dEFT to the experimental data [HH - H2] . 

be improved. We try to incorporate the formulae of dEFT. 

Figure shows the total cross-section calculated by using G4dEFT (red solid line), to¬ 
gether with the experimental data [54fl42] . and the evaluated nuclear data from ENDF/B- 
VII.1 |l3j and JENDL/PD-2004 |11]. ENDF/B-VII.l does not reproduce some experimental 
data around the peak, but our G4dEFT result shows good agreement with data not only 
around the peak, but also over the energy range considered. 

The ratios of the total cross sections obtained by using G4dEFT (denoted by G) to those 
from experiments (denoted by E) [MH12] are shown Figure The results from G4dEFT 
are consistent with the experimental data for the photon energy below 15 MeV within 
experimental errors. For the region 15 < Ej < 30 MeV, the ratios are largely located in the 
range from 1.0 to 1.2. There are large discrepancies (factors of ~ 2) for E^ = 25.4 MeV and 
E-y = 25.7 MeV, but experimental errors are also large. For the photon energy E^ above 30 
MeV, four experimental data are shown and agree with those from G4dEFT within about 
20% error. We note that, in principle, the pionless theory will break down for momenta 
larger than the mass of pion. In this respect, it may make sense to compare the result with 
data for the photon energies below about 20 MeV. One can see that the calculated results 
agree well with the data in this region. 


4 










FIG. 3: (Color online) Differential cross-section (j{9) as a function of cos 9 multiplied by a factor 
100 /(Tt, where at denotes the total cross-section of the unpolarized photons at the energies (a) 

= 19.8 MeV and (b) = 38.6 MeV. The open squares represent the experimental data taken 

from Ref. 46, and the solid lines denote the results obtained from G4dEFT. 

B. Differential Cross-Sections 

In order to investigate the qnality of the model in more detail, we compare the differential 
cross-sections obtained from G4dEFT with the experimental data available in literature or 
the EXFOR database |1S]. In Fig. we present the differential cross-sections at photon 
energies (a) = 19.8 MeV and (b) = 38.6 MeV as functions of the polar angle 6 in the 

c.m. frame. Open squares are the experimental data from Ref. 46, and solid lines represent 
the results of G4dEFT. For E^ = 19.8 MeV, G4dEFT agrees with the experiment data 
within the error bars, but for E^ = 38.6 MeV, there are systematic underestimates in the 
backward angles. 

Figure 1^ shows the differential cross-sections at 6 *iab = 45°, 75°, 90° and 105° as functions 
of the incident photon energy E^, where 6 'iab is the polar angle in the laboratory frame. 
Open squares and circles represent the experimental data taken from Ref. 47 and Ref. 48, 
respectively. Solid lines denote the results obtained from G4dEFT. In the energy range 
E^ < 30 MeV, the results of G4dEFT are mostly within the error bards regardless of the 
detection angles. 

Some experimental works report the ratio of differential cross sections by taking ^lab = 90° 
as a reference angle [IHIEO]. These works provide data sets independent of the differential 
cross-sections in Figs. |^and|^ In Fig. we show differential cross-sections at 6 *iab = 45°, 
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FIG. 4: (Color online) Differential cross-sections at 6\ah = 45°, 75°, 90° and 105°. Open squares 
and circles correspond to the experimental data in Ref. 47 and Ref. 48, respectively, and the values 
are taken from the EXFOR database |45] . Solid lines are the results obtained from G4dEFT. 

135°, and 155° divided by those at 6*iab = 90°. Experimental data are denoted with symbols, 
and G4dEFT resnlts are depicted in solid lines. On the average, onr resnlts deviate from the 
data of Stephenson et al., Ref. 49 by abont 22%. Comparing the G4dEFT resnlts to the data 
by Birenbanm et al., Ref. 50, we have abont 6% discrepancy in the energy range 7 < E-y < 9 
MeV for ^lab = 45°, bnt similar disagreement as that with Stephenson et al. appears at other 
angles and energies. G4dEFT cross-sections underestimate the experimental data from Ref. 
50 by about 15% and 31% for 6*iab = 135° and ^lab = 155° on average, respectively. On the 
other hand, a comparison with the data of Sawatzky ISB gives relatively good agreement 
at energies E.y < 6 MeV. Comparison with the data of Blackston isa shows agreement at 
Oiah = 45° with about 4.3% error, but the disagreements at ^lab = 135° and 155° are not 
resolved. 
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FIG. 5: (Color online) The symbols denote the experimental data in Refs. 49-52, and the actual 


numbers are taken from the EXFOR database |45| . 


C. Neutron Yields from CD2 Photodisintegration 

So far, we have focused on checking the applicability of G4dEFT by comparing our 
results for dj —?• np to experimental data. As an application of G4dEFT, we attempt to 
predict the neutron yields from the 7 -|- GD 2 process, which has not yet been explored either 
experimentally or theoretically. In this process, neutrons can be emitted from both carbon 
and deuteron. As shown in the dy —)■ np results in Fig. [T| different results can be anticipated 
from G4 and G4dEFT. One can think of gaseous (D 2 ) or liquid (D 2 O) targets, but taking 
into account the advantages of the solid target, such as easy handling and large reaction 
rates, we chose GD 2 as a simulation material. 

A schematic geometrical setup of our simulation is depicted in Fig. We assume a 
pencil beam for the incident photons. The GD 2 target with a density of 1 g/cm^ is modeled 
as a cylinder 1 cm in diameter and 1 cm in thickness. The scoring region with the shape of 
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FIG. 6: (Color online) Schematic diagram of the simulation setup. 

a spherical shell surrounds the target. The inner and the outer radii of the scoring region 
are chosen to be 100 cm and 100.1 cm, respectively. The target area is in vacuum, and the 
thickness of the scoring region is chosen as 0.1 cm arbitrarily for convenience and does not 
affect the results. 

For the EM processes and hadronic interactions of hadrons in our simulations, we use 
GEANT4 physics lists, “G4EmLivermorePhysics” and “G4HadronPhysicsQGSP_BIG_HP” 
classes, respectively.^ As mentioned in Sec. II, GEANT4 uses G4EmExtraPhysics for pho- 
tonuclear interactions. G4EmExtraPhysics can be applied to carbon materials, but as we 
have shown, it is not valid for the deuteron in the low-energy region. To see the effect of the 
deuteron, we simulate the 7 -|- GD 2 process by considering two cases. One is the case using 
only G4EmExtraPhysics (G4), and the other is the case using both G4EmExtraPhysics and 
G4dEFT (G4 + G4dEFT). 

Figure shows the distribution of neutrons scored in the shell as a function of the neutron 
energy En for (a) = 9.8 MeV, and (b) E^ = 19.8 MeV. In each plot, the results of G4 

(G4-|-G4dEFT) are denoted as hlled circles (open squares). Let us analyze the results of G4 
first. As shown in Fig. the cross section of dy —?• np is null in G4, so the neutrons for 
G4 in Fig. are emitted only from carbon nuclei in GD 2 . For E^ = 9.8 MeV, very small 

^ There are many different pre-defined physics constructors in GEANT4. To check the dependence on the 
constructors, we simulated both the photon flux in the CD 2 target and the neutron yields generated by 
c ?7 —>■ np reaction with different GEANT4 physics constructors, separately, but the differences due to the 
constructors were almost negligible because the simulated observables discussed in this work came from 
the dy —)■ np reaction by using G4dEFT. 
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FIG. 7: (Color online) Expected energy spectra of neutron yields due to photons of energies (a) 
9.8 MeV and (b) 19.8 MeV incident on a CD 2 target. The black solid line with filled circles and 
the red dotted line with open squares represent the results obtained by using GEANT4 with G4 
and G4 + G4dEET, respectively. 

numbers of neutrons are emitted up to ~ 1 MeV, but more energetic neutron are not 
produced. In the target, we assume natural carbon, in which about 99% is and 
the remaining 1% is mostly The threshold for the ^^ 0 ( 7 ,n) reaction is 18.737 MeV, 
so the black circles in Fig. (a) denote the neutrons emitted from only. Because the 
fraction of is very small in nature, the production of neutrons is also very small. Due 
to energy conservation, the maximum energy of the emitted neutrons is limited. This is the 
reason the neutron production is truncated at ~ 1 MeV. For = 19.8 MeV, the photon 
energy exceeds the threshold of the ^^C( 7 ,n) reaction, so we have a substantial number of 
neutrons in the energy range = 0 ~ 1 MeV. The number of neutrons drops abruptly 
around ~ 1 MeV, which is also due to the conservation of the energy of the neutrons 
emitted from The black circles indicate the number of neutrons emitted from at 
energies above ~ 1 MeV. Here again, the small numbers in the range En > 1 MeV reflect 
the small fraction of in nature. 

Now, let us discuss the results of G4+G4dEFT. In this model, the neutrons are con¬ 
tributed by both carbon and deuteron. The neutron yield in the G4-|-G4dEFT model is 
similar to that in the G4 model in the ranges < 1 MeV for E^ = 9.8 MeV and E^ < 7.4 
MeV for E^ = 19.8 MeV. However, signihcant differences occur around En ~ 3.7 MeV for 
E-^ = 9.8 MeV and around En ~ 8.7 MeV for E^ = 19.8 MeV, where high and broad peaks 
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FIG. 8: (Color online) Expected angular distributions of neutron yields in the G4+G4dEFT model. 
The black dotted line with squares and the red solid line with circles represent the results for 9.8 
MeV and 19.8 MeV photons, respectively. 

appear. The number of neutrons produced by the d'y —?• np reaction is obtained by making 
use of the “GdUserSteppingAction” class. The ratio of the neutron number for = 9.8 
MeV to that for = 19.8 MeV is obtained as 2.293 ± 0.046. This number is close to the 
ratio of the total cross-sections of the dj —)■ np process, 2.287, for E^ = 9.8 MeV (1.44 mb) 
and for E^ = 19.8 MeV (0.63 mb). The similarity of two ratios is proof that the peaks in 
Fig. [^originate from interactions of photons and deuterons. 

In Fig. we show the angular distribution of the neutron yield. The maximum values 
are located around ^lab = 90°, which is in agreement with the behavior of the differential 
cross-section of the d'y —)■ np process, Fig. This coincidence proves again that most of 
the neutrons in the peaks are emitted from the deuterons in CD 2 . The angular distribution 
for E^ = 19.8 MeV is relatively broader than that for E^ = 9.8 MeV, which is one of the 
reasons we have a peak broader for E^ = 19.8 MeV than for E^ = 9.8 MeV in Fig. 

The distribution of the emitted neutrons in terms of the energy at 6*iab = 90° is plotted in 
Fig. 1^ Very sharp peaks are seen at En ~ 3.7 MeV for 77..,, = 9.8 MeV and at 77^ ~ 8.7 MeV 
for E^ = 19.8 MeV. If one treats the scattering process d'y —>■ np classically and calculates 
the energy of outgoing neutrons in the laboratory frame, one can easily obtain E^ = 3.7 
MeV and 8.7 MeV for 77..^ = 9.8 MeV and 19.8 MeV, respectively. The coincidence of these 
values and the positions of the peaks signal the dominance of dy reactions in the emission 
of neutrons. 
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FIG. 9: (Color online) Expected neutron energy spectra at 0iab = 90° produced by photons on a 
CD 2 target in the G4+G4dEFT model. 

IV. SUMMARY 

With the development of effective held theories in nuclear physics, the accuracy in de¬ 
scribing and understanding few-nucleon systems is being improved unprecedentedly. The 
improvement will make it possible to have more precise results for simulations in which 
few-nucleon processes will contribute non-negligibly. Our work started from the question of 
how we could embed the progress of nuclear theory for few-nucleon systems in the simula¬ 
tions of nuclear processes. As a hrst step, we checked the validity of G4EmExtraPhysics, 
which accounts for the photonuclear interactions in GEANT4 (vlO.l) by calculating the 
total cross-section of the dy —)■ np process. We found that the G4EmExtraPhysics model 
gave null result for the total cross-section of the process at energies below the pion thresh¬ 
old. In order to resolve the problem and enhance the predictive power of simulations, we 
incorporated analytic results of a dibaryon effective held theory into GEANT4. 

The new model has been checked thoroughly by calculating the total and the differential 
cross-sections of the dy —)■ np process for incident photon energies up to 70 MeV and by 
comparing the results with available experimental data. G4dEFT gives results within the 
experimental error bars for below 15 MeV. For photon energies above 15 MeV, the results 
from G4dEFT show reasonable agreement with those from experiments within 20%, even 
though the energy region is out of the valid range for the theory. 

The model has been applied to calculate the yield of neutrons emitted from GD 2 bom¬ 
barded by photons. Because no experimental datum is available for the observable, only a 
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comparison between the new and the existing models was made. The two models show a 
critical difference in the prediction of the neutron yield. The measurement, if performed, 
will provide criteria about the accuracy and the extent to which the model is valid and 
applicable. 
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